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RESEARCH  OF  THE  PROJECT  FOR  ADJUSTING  THE  THRUST  OF  A 
SOLID  ROCKET  PROPULSOR 


Xu  Wengan 


ABSTRACT 

This  paper  presented  a  proposed  plan  for  a  solid 
rocket  motor  with  adjustable  propellant.  It  was  real¬ 
ized  through  the  use  of  a  negative  pressure  exponent  pro¬ 
pellant.  In  this  paper,  the  feasibility  and  characteris¬ 
tics  of  the  adjustment  of  this  technique  were  discussed. 


It  also  pointed  out  the  other 
pressure  exponent  propellant. 
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I .  INTRODUCTION 


In  the  area  of  thrust  control  of  a  solid  rocket  engine,  a  large 
amount  of  work  has  been  performed  both  here  and  abroad.  As  a 
result,  directional  control  systems  and  thrust  termination  mechan¬ 
ism  of  various  principles  such  as  the  "deflected  flow  ring",  "rotat¬ 
ing  nozzle",  "vibrational  nozzle",  "soft  nozzle",  "liquid  suspension 
nozzle"  and  the  "secondary  nozzle"  for  liquids  or  combustible  gases 
were  developed.  However,  with  regard  to  the  adjustment  of  the  mag¬ 
nitude  of  the  thrust,  there  is  not  yet  a  good  feasible  method. 

Although  some  people  had  proposed  some  possible  means  such  as  the 
"layer  engine"  and  the  injection  of  a  certain  liquid  to  the  combus¬ 
tion  chamber  and  some  preliminary  attempts  [1,2]  were  made  using  a 
propellant  with  a  positive  pressure  exponent  n,  which  is  greater  than 
zero  and  close  to  1,  and  a  mechanical  rotating  valve  or  an  adjustment 
cone  in  the  nozzle,  however,  due  to  the  complicated  structure  and  the 
poor  weight  and  performance  characteristics,  these  methods  do  not  have 
an  optimistic  future. 

This  paper  is  an  attempt  to  use  a  new  propellant  with  a  less  than 
zero  pressure  exponent  whose  absolute  value  is  greater  than  one.  For 
ease  of  discussion,  the  author  named  it  the  "high  negative  pressure 
exponent"  propellant  as  the  concept  to  solve  this  problem. 

II.  THE  FEASIBILITY  OF  USING  HIGH  NEGATIVE  PRESSURE  EXPONENT 
PROPELLANT  TO  REALIZE  THE  ADJUSTMENT  OF  THE  MAGNITUDE  OF 
THRUST  AND  ITS  ADVANTAGES 

Why  do  we  have  -to  search  for  this  type  of  propellant?  Can  this 
propellant  provide  the  adjustment  of  thrust?  What  are  its  advantages? 
Obviously,  these  are  the  questions  to  be  answered  first.  As  it  is 
commonly  known,  when  the  combustion  speed  follows  the  r  =  op,*  law,  a 
sonic  flow  always  exists  at  the  throat  of  the  nozzle.  The  equilibrium 
pressure  and  the  corresponding  propulsion  in  the  combustion  chamber  of 
a  solid  rocket  engine  can  be  expressed  by  the  two  following  formulas: 

j _ 

Pc  =  (  (1) 
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From  the  above  equation,  in  order  to  realize  the  adjustment  of 
thrust,  the  only  effective  way  is  to  change  the  two  parameters  A^ 
and  .  Among  them,  the  only  parameter  to  realize  the  random  change 
using  a  simple  method  (i.e.,  the  type  of  parameter  which  does  not 
follow  a  pre-determined  program  but  changes  according  to  the  need 
according  to  the  circumstance)  is  mainly  the  throat  area  of  the  nozzle 
At.  Therefore,  we  will  use  the  following  equation  obtained  from  the 
simultaneous  solution  to  (1)  and  (2)  as  the  basic  equations  to  be  dis¬ 
cussed  in  this  paper: 

Fe-K,A,  n~1 
F  \  »- 


Mfc) 

PC  =  K,FX' ' 


(3) 

(4) 

(5) 


where 


K  i—(P  r'C*’Ab-a') 
K,  =  Cr-Kt 


1 

i — 1» 


K,=  K,  •  K,’ 


when  the  propellant  is  chosen  and  the  combustion  area  Ab  is  a  constant 

obviously  it  is  allowable  to  use  K.  as  a  constant  to  analyze  the 

1  A 

problem.  Furthermore,  C„  is  also  a  function  of  _b  and  P  .  However, 

r  A  c 

within  the  realm  of  this  discussion  it  can  be  t  treated  as  a  con¬ 
stant  approximately.  Therefore,  in  the  following  discussion,  K2  and 
are  also  treated  as  constants  in  approximation.  Thus,  by  taking 
the  logarithmic  differentiation  on  both  sides  of  the  above  equation, 
we  get 


From  equations  (5)  or  (8),  we  know  that  the  adjustment  of  the 
magnitude  of  thrust  will  cause  a  pressure  variation  inside  the  com¬ 
bustion  chamber.  However,  when  the  pressure  variation  becomes  totally 
unacceptable  to  the  practice,  the  adjustment  of  the  magnitude  of  the 
propulsion  becomes  an  empty  thought.  Therefore,  it  is  necessary  to 
minimize  the  variation  of  combustion  chamber  pressure  during  the  adjust 
ment  of  propulsion  (consequently,  the  structural  strength  of  the 
engine  can  be  sufficiently  utilized).  At  this  time,  the  absolute  value 
of  n  should  be  as  large  as  possible.  For  example,  when  the  pressure 
exponent  n  =  0.75,  in  order  to  increase  the  thrust  by  eight  or  four 
times,  the  pressure  must  be  increased  by  16  or  6.3^  times.  These 
changes  are  so  huge  that  the  weight  performance  of  the  engine  is 
greatly  reduced  in  order  to  ensure  the  engine  to  function  under  the 
maximum  pressure.  This  is  the  fatal  weakness  of  the  experimental  and 
development  plans  reported  abroad.  However,  when  n  =-2,  to  increase 
the  thrust  by  the  same  factors,  the  pressure  is  reduced  to  about  1/3.83 
or  1/2  of  the  original  pressure.  Obviously,  the  amplitude  of  pressure 
variation  inside  the  combustion  chamber  is  significantly  reduced. 

On  the  basis  of  the  above  consideration,  obviously,  the  large  |n| 
is  the  better.  However,  in  order  to  let  the  engine  operate  stably, 
the  n  value  with  an  absolute  value  greater  than  1  cannot  be  positive. 
This  is  because  when  the  flow  rate  of  the  engine  and  the  formation 
rate  of  the  combustible  gas  are  equal,  it  will  work  under  an  equili¬ 
brium  pressure  (pcb).  This  equilibrium,  just  as  any  other  equilibrium 
problem,  also  has  the  stable  equilibrium,  unstable  equilibrium  and 
occasional  equilibrium  conditions.  In  other  words,  if  some  factor 
happens  to  affect  the  pressure  of  the  combustion  chamber  to  deviate 
from  the  equilibrium  ?cb,  then  is  there  a  tendency  to  decrease  or 
increase  this  deviation?  Or  is  it  reaching  equilibrium  under  new 
conditions?  These  are  the  differences  of  the  three  types  of  equilibria. 
For  an  engine  to  work  stably,  obviously  we  must  have  a  stable  equili¬ 
brium.  The  mathematical  expression  of  this  analysis  may  be: 


where 


therefore 


because  at  the  equilibrium  (Pf  —  Pc)AtaP"  = 
point  we  have 

Ry  substituting  into  the  f-  "  )  =  — -p^ —  Ln  -  ]  ] 

above  equation,  we  get  '  '  "  ' 


also  because 


rir *'•  '•  4 

K.  C*  ^ 


Therefore,  we  have  to  let-— ^ — /  JL£s—)<q 

dp,  V  / 


(Note  1)  "ip- (-*&-)<>  is  the  common  mathematical  expression  _*LL~>o. 

r 


because  the  former  can  be  converted  to 


dp,  ( 


It  is  necessary  to  let 
which  is 


.  which 

)<°  ^  T~E-y  <° 
V  dPr  ) 

+•.  >o 

dp\ 

n  -  1  <  0 
n  <  1 


This  can  be  intuitively  observed  from  Figure  1. 

For  this  reason,  in  order  to  minimize  the  combustion  chamber 
pressure  variation  during  the  adjustment  of  the  magnitude  of  the 
thrust,  we  should  use  a  "high  negative  pressure  exponent”  propellant. 


Simultaneously,  from  equations  (1)  and  (2)  as  well  as  the  logar¬ 
ithmic  differentiation  of  both  sides  of  equation  (1),  we  get 

dp.  If  dP,  dC*  .  dA>  dA,  d.  "I  ,  . 

~p,  ‘  ~  T-n  LV  C’  X  ~A, "  ~a  J 

We  can  see  that  the  fluctuations  of  parameters  such  as  p,,C*.A.^*.° 
caused  by  a  certain  factor  will  cause  the  fluctuations  of  the  combustion 
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Figure  1  Figure  2 

3 — formation;  4 — flow  out;  7 — flow  out;  8 — flow  out; 

5 — formation;  6 — formation  9 — formation;  10 — formation 

chamber  pressure  p  and  the  propulsion  force  F.  However,  when  the 

V 

relative  fluctuation  of  the  former  is  fixed,  the  relative  fluctuation 
of  the  latter  decreases  with  decreasing  n  value.  In  practice,  the 
maximum  fluctuation  of  the  above  parameters  (such  as  Pp.C*%Ai%A,.a  ,  etc.) 
is  the  combustion  speed  coefficient  a.  However,  from  the  actual 
charge  production  requirements,  it  is  frequently  required  to  fluctuate 
in  a  wider  range.  For  example,  a  heterogeneous  propellant  has  two 

A  a 

typical  formulations  and  their  —  are  regulated  to  be  about  0.07  and 

cL 

0.02  respectively.  This  is  already  very  difficult  technologically. 

However,  if  n  =  0.75>  together  with  the  fluctuations  of  other  parameters 

has  already  reached  as  high  as  0.28  and  0.08,  respectively.  If 
rc  APc 

n  =  -2,  then  — p^-  can  be  reduced  to  0.023  and  0.007,  respectively.  Its 
advantage  is  obvious.  It  can  indeed  significantly  reduce  the  fluctua¬ 
tion  of  combustion  chamber  pressure  and  the  distribution  of  the  propul¬ 
sion  of  the  engine.  In  addition,  the  effect  of  on  is 

similarly  reduced  significantly.  Therefore,  it  allows  the  use  of  a 


high  surface  ratio  charge  design  to  improve  the  filling  density  of 
the  engine. 

8 

Summarizing  the  above,  as  shown  in  Figure  2,  after  using  the  new 
"high  negative  pressure  exponent"  propellant,  it  is  possible  to  adjust 
the  thrust  within  a  relatively  large  amplitude  under  the  condition 
that  the  present  structure  of  the  present  solid  rocket  engine  might 
not  have  to  be  changed.  In  addition,  it  has  the  advantages  of  low 
combustion  chamber  pressure  flutuation,  small  engine  propulsion  dis¬ 
tribution  and  high  filling  density. 

Obviously,  in  order  to  turn  the  adjustment  of  the  magnitude  of  the 
thrust  into  reality,  we  must  allow  the  throat  area  to  be  adjustable  by 
a  large  amplitude.  This  can  be  seen  from  equations  (4)  and  (7).  We 
can  see  that  when  n  <  0  and  |n|  is  greater  than  1,  in  order  to  obtain 
the  needed  adjustment,  we  must  be  able  to  adjust  the  nozzle  throat  area 
significantly.  As  shown  in  the  two  examples  given  above,  for  an 
increase  in  the  propulsion  force  of  eight  times,  At  must  be  reduced 
by  a  half  for  n  *  0.75*  However,  for  n  *-2,  At  must  be  increased  to 
22.6  times  of  the  original  value.  This  is  not  possible  for  the  usual 
mechanical  adjustment  of  the  throat  area  of  the  nozzle.  However,  if 
a  flow  controlled  "vortex  value"  is  used,  there  it  is  no  problem  at 
all.  Furthermore,  with  the  increasing  energy  of  the  propellant,  it  is 
necessary  to  rely  on  the  flow  control  technique  to  adjust  the  throat 
area  of  the  nozzle  (i.e.,  the  flow  rate)  when  the  nozzle  throat  is 
working  in  such  hostile  conditions.  Therefore,  in  theory,  this 
requirement  is  not  an  obstacle  to  the  realization  of  thrust  adjustment 
based  on  this  method.  On  the  contrary,  it  provided  the  feasibility 
for  the  "microadjustment"  of  the  thrust.  Therefore,  it  can  be  pre¬ 
dicted  that  it  will  have  the  ideal  static  characteristics  and  accuracy 
as  a  link  of  the  automatic  adjustment  system. 

III.  CAN  THE  "HIGH  NEGATIVE  PRESSURE  EXPONENT"  PROPELLANT  BE  OBTAINED? 

Due  to  w.,c  1?  ited  work  available,  the  author  cannot  perform  a 
quantitative  analysis  based  on  the  combustion  mechanism.  Now,  the 
following  points  are  presented  as  the  basis  of  discussion: 
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1.  Not  only  in  theory,  but  also  in  practice  there  are  n  <  0 
propellants  in  existence.  For  example,  in  a  certain  formulation, 
within  the  pressure  range  from  50-70  kg/cm  .  the  conditon  that  the 
pressure  exponent  n  =-0.825  exists. 

In  the  reports  abroad,  the  so-called  "mesa"  effect  exists  in  the 
double  base  platform  propellanc  and  heterogeneous  propellant  with 
n  <  0  [3,5].  It  was  further  reported  that  heterogeneous  propellants 
with  n  <  0  or  even  n  =-4.0  existed  [4].  This  shows  that  it  is  feas¬ 
ible  and  practical  to  obtain  a  "high  negative  pressure  exponent"  pro¬ 
pellant  . 

2.  In  the  foreign  literature  [1,4, 5, 6],  the  mechanism  of  the 
"mesa  effect"  regarding  the  combustion  of  heterogeneous  solid  pro¬ 
pellants  have  been  described. 

(1)  From  the  micrographs  of  the  combustion  surface  after  the 
fire  is  ceased  in  experiments,  we  can  see  that  under  low  pressures 
(less  than  30  kg/cm  )  the  ammonium  perchlorate  (Ap)  crystals  are 
protruding  on  the  surface.  This  is  because  under  such  conditions  the 
thermal  decomposition  of  the  binder  is  faster  than  that  of  Ap. 

While  under  high  pressures  (greater  than  30  kg/cm  ),  the  results  are 
the  reverse.  There  are  holes  indicating  that  the  Ap  crystals  are  con¬ 
sumed  very  rapidly  [1]. 

(2)  The  pU  binder  is  different  from  binders  such  as  CTPB.  It 

is  already  very  fluid  before  reaching  the  rapid  decomposition  tempera¬ 
ture.  However,  the  latter  is  a  high  viscosity  fluid  with  bubbles. 
Under  .high  pressures,  a  local  cease  fire  effect  will  appear  [1,4, 5, 6] 
when  the  highly  fluid  binder  flows  into  and  covers  the  holes  of  the  Ap 
crystals  on  the  combustion  surface. 

(3)  It  was  observed  using  a  movie  camera  that  there  was  the 
irregular  pulsing  combustion  effect  with  local  cease  fire  in  the  flat 
form  propellants.  The  magnitude  of  the  fraction  of  local  cease  fire 
with  respect  to  the  total  area  determines  the  extent  or  average  com¬ 
bustion  speed  decrease  on  the  full  combustion  plane.  However,  the 


dependence  of  this  fraction  on  pressure  will  cause  the  reduction 
of  the  pressure  exponent  n,  which  initiates  the  "mesa  effect"  [1,4,51 
with  n  <  0 . 


(4)  With  the  addition  of  an  easily  fuseable  salt  such  as 
(NH^)2SOi|  as  the  combustion  speed  adjuster,  it  can  cause  a  similar 
local  cease  fire  affect  or  create  a  deep  thermal  potential  base  and 
the  suppression  effect  of  NH^  on  the  decomposition  of  NH^CIC^  [4]. 

3.  Reference  [4]  presented  a  physical  model  of  the  combustion 
process  of  this  type  of  propellant  which  was  partially  verified  by 
experimental  observations.  An  experiment  of  the  engine  using  80 
pounds  of  propellants  with  n  =  -2.5  was  also  reported. 

IV.  PRELIMINARY  ANALYSIS  OF  SOME  ADJUSTING  CHARACTERISTICS 


For  the  large  circuit  of  the  automatic  propulsion  adjustment  sys 
tern,  the  adjustment  of  the  engine  is  only  one  of  the  links.  Its 
adjusting  quality  obviously  must  be  analyzed. 


For  convenience,  let  us  assume  that  the  combustion  speed  still 
follows  the  ■  =arr  law  during  the  adjustment  process  that:  there 
is  a  sonic  flow  at  the  nozzle  throat  at  all  the  times,  the  pressure 
inside  the  combustion  chamber  is  equal  everywhere  and  it  is  equal  to 
the  nozzle  inlet  pressure.  Thus,  from  the  continuity  equation  of 
the  combustable  gas: 


dP, 

it 


r:r*: 


I  (Pt  -i\)A.up'~ 


C‘ 


(10) 


From  the  logarithmic  differentiation  with  respect  to  equation 

(2)  and  by  assuming  Cp  is  a  constant,  we  get 

dF  _  dp,  _  dA, 

"  F  pr  A, 

or 

1  dF  1  / -i -1  'I 

F~  it  ~  V  '  it  A,  tit  UiJ 

By  solving  equations  (2),  (10)  and  (11)  simultaneously,  we  can 
obtain  the  dynamic  function  of  this  link  as: 


3  dF  i-  /  F  \n-'  i  .  , 

>-*• -*••<•** 


where 


r*r*! 

A.=  .  ■—  (np-rc)A,aCy 

»  . 


When  the  input  is  a  unity  step  function,  i.e.,  )  (A^  is  the 

throat  area  before  the  unity  step  function),  obviously,  when  t  >  0, 
u(i)=i,  dA'  —  q  .  At  this  time,  the  differential  equation  (12)  is: 

dt 


where 


d?-  +  K\F  =  K\F' 

dt 


l’  •  _ _ A.  , 


The  obtained  unity  transition  differential  equation  (13)  is  the 
classical  Bernoulli  equation.  Hence,  through  a  variable  separation 
method,  we  get:  when  n  ^  0  _i  ,  ,  ,..r-  ,  r 

t  =  — —  ln(  A A  ,r  )  + 

(,  1  —  n)K  i 


From 


the  initial  condition,  t  =  0+F=F(0),  we  get  C,  =  (Y-ny/r" in  K',-K‘2IF(Q)2 


(1—  n)Kj 


.  A';-/ficF(o)i: 

n  KV-KiF " 


Therefore 


tr.  f(f)=;  {l  — A,exPL(w— 

Where  , 

*■-(-*!  )*=■' 

Also  because  (n-1)  <  0,  therefore,  from  equation  (14)  we  can 

r  p/  \ 

get  when  i-+  oo  ,  F(og)  =  ;  .  Simultaneously,  we  can  get  /.;  =  ilir: '  - 

L  •  ( 0 )  j 

Therefore,  equation  (14)  can  be  rewritten  as: 


F(O-F(0)Jl-[l-(-^-)),’"]«Pt(n-l)A'i<]'»  «-■ 

1.  In  the  transition  period:  from  the  definition 

!F(r,.„)— F(°°)J_  — Q  05  can  get 
F(  oo) 


(15) 


10 


To-os — ;  m.'r 

(n-l)A, 


’“l1  ±0.05)" 
V  F( 0)  j'-" 

-F(°o)j 


where  the  +  and  -  signs  correspond  to  the  0<«<]  and  n<o  conditions, 
respectively . 


When 


the  relative  adjustment  amount  -^7 -{^-—-*100#  =  10% 


F(  oc  ) 


corresponding  to  the  conditions  of  n  =  0.75  and  n  =  -2,  tn  nt-  are 

143  0213  U*Uj 

about  —g.  —  and  ,  respectively.  This  also  means  that  the 

transition  period  of  the  latter  is  about  1/7  of  the  former. 

2.  The  over  adjustment  quantity:  From  equation  (14),  we  can 

get  Fmax  =  Therefore,  according  to  the  definition,  the  over 

adjustment  0  =  F  _ ioo#=0. 

F(oo) 

3.  Static  difference:  From  the  given  F#  =  /\.*Cf*M,,+ii(0] 

and  p.,-k->,)-c-74i#  • we  can  Eet 

L  .  Therefore,  according  to  the  defini¬ 

tion,  the  static  differencet(oo)~l^^^”\ioof> «o(the  P  above  is  the 

M*  g 

static  equilibrium  pressure  corresponding  to  a  given  propulsion). 


Therefore,  when  the  pressure  changes  rapidly,  the  combustion 
speed  is  no  longer  a  function  of  pressure  alone.  It  is  also  a  func¬ 
tion  of  the  pressure  variation  rate  such  as  r— o/>;  (l  +  ~£~) 

(assuming  that  this  empirical  formula  is  applicable  to  the  n  <  0 
condition).  Similarly,  the  dynamic  equation  corresponding  to  a  unity 
step  process  can  be  obtained  as: 

dF  _  KrF~"-KtF''' 

~d,  K>F*4'  —  KT 

when  **1. 0  ,  its  solution  is 

dF  \ 

t  =  ln(  K\ -  K\F 1  *■)  -  I  XtF*^r  ) 

(1  —n)Kt  ■  * 

In  the  meantime,  there  are  experiments  which  showed  that  the 
combustion  pressure  exponent  n  of  a  pressure  transient  process  is 

9. 

greater  than  the  pressure  exponent  n  under  a  gradual  change  [7]. 
Therefore,  the  assumption  used  above  based  on  the  law  in  the 

adjustment  process  is  obviously  different  from  reality.  This 
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difference  can  be  observed  from  the  experimental  regularity  shown 
in  Figure  (2)  in  [**].  Therefore,  the  above  analysis  is  only 
preliminary. 
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